Observation of a Turbulence-Induced Large Scale Magnetic Field 
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An axisymmetric magnetic field is applied to a spherical, turbulent flow of liquid sodium. An 
induced magnetic dipole moment is measured which cannot be generated by the interaction of the 
axisymmetric mean flow with the applied fleld, indicating the presence of a turbulent electromotive 
force. It is shown that the induced dipole moment should vanish for any axisymmetric laminar flow. 
Also observed is the production of toroidal magnetic field from applied poloidal magnetic field (the 
oj-effect). Its potential role in the production of the induced dipole is discussed. 
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Many stars and planets generate their own nearly- 
axisymmetric magnetic fields. Understanding the mech- 
anism by which these fields are generated is a problem 
of fundamental importance to astrophysics. These dy- 
namos are sometimes modeled using two components: 
a process which generates toroidal magnetic field from 
poloidal field and a feedback mechanism which reinforces 
the poloidal field The first process is easily modeled 
in an axisymmetric system: toroidal differential rota- 
tion of a highly-conducting fluid sweeps the pre-existing 
poloidal field in the toroidal direction creating toroidal 
field. This phenomenon, known as the w-effect, is effi- 
cient at producing magnetic field and has been observed 
experimentally |^ Ql . The second ingredient to the 
model is more subtle, as toroidal currents must be gen- 
erated to reinforce the original axisymmetric poloidal 
field. Cowhng's theorem [3 excludes the possibility of 
an axisymmetric flow generating such currents so some 
symmetry-breaking mechanism is required. 

The usual mechanism invoked 0| is a turbulent electro- 
motive force (EMF), £ — (v xh^, whereby small scale 

fluctuations in the velocity and magnetic fields break 
the symmetry and interact coherently to generate the 
large scale magnetic field. This EMF is sometimes ex- 
panded in terms of transport coefficients about the 
mean magnetic field: £ — aB -I- /3V XB-I-7XB; a is 
characterized by helicity in the turbulence, (3 by en- 
hanced diffusion and 7 by a gradient in the intensity 
of the turbulence, a is of particular interest as it re- 
sults in current fiowing parallel to a magnetic field, and 
when coupled with the w-effect can generate the toroidal 
currents needed to reinforce the poloidal field. 

Experimental evidence for mean-field EMFs (such as 
the a-effect) in turbulent fiows has been scarce. Three 
experiments, relying on a laminar a-effect, have gener- 
ated an EMF Q and dynamo action 0, , but heavily- 
constrained flow geometries were used to produce the 
needed helicity; the role of turbulence was ambiguous. 
Experiments with unconstrained flows have provided ev- 
idence for turbulent EMFs, though not the turbulent a- 




FIG. 1: Schematic of the Madison Dynamo Experiment show- 
ing a cut-away view of the sphere, impellers, external field 
coils, surface and internal Hall probes. 

effect. Rcighard and Brown 11] have attributed a mea- 
sured reduction in the conductivity of a turbulent flow of 
sodium to the /3-effect. Petrelis et al. have observed 
distortion of a magnetic field similar to an a-effect (cur- 
rents generated in the direction of an applied magnetic 
field) and postulate that turbulence may be responsible 
for disagreement between a laminar model and observa- 
tions. Not all liquid-metal experiments have had such 
results: Frick et al. have reported that the mean 
fiow accounts for all magnetic fields in their torus-shaped 
gallium experiment, and Peffley, Cawthorne and Lath- 
rop Q have observed no such effects. It should also be 
noted that an a-effect has been observed in the core of 
magnetically-confined plasmas 0, 0| . 

In this Letter we report measurements of the magnetic 
field induced by applying an axisymmetric magnetic field 
to a turbulent, axisymmetric fiow of liquid sodium. An 
induced dipole moment is measured which cannot be gen- 
erated by the mean flow, indicating the presence of a 
turbulent EMF. 

The study is conducted in the Madison Dynamo Ex- 
periment, aim diameter stainless steel sphere containing 
liquid sodium. As shown in Fig.^ two drive shafts enter 
the sphere through each pole and drive 30.5 cm diame- 
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FIG. 2: Upper half: color contours of induced toroidal mag- 
netic field, B^{s, z), measured by sets of internal Hall probes, 
for Rmtip = 100. Induced poloidal flux surfaces, ^(s,z), 
are in black. The positions of the internal Hall probes are 
indicated with dots. The cylindrical axis of symmetry is hor- 
izontal. Lower half: velocity field measured in a water model 
of the Madison Dynamo Experiment, for an impeller rota- 
tion rate of 16.7Hz. Contours of toroidal flow, v^{s,z), are 
in color and poloidal stream function, <l>(s,z), are in black. 
The arrows indicate the direction of the poloidal flow, and 
the rectangles indicate the positions and size of the impellers 
which drive the flow. 



ter impellers which generate an axisymmetric mean flow. 
The shafts are coupled to two 75 kW motors which are 
independently controlled by variable-frequency drives. 
The radial component of the magnetic field is measured 
by an array of 74 temperature-compensated Hall probes 
mounted to the sphere's surface, allowing resolution of 
spherical harmonic components of the external magnetic 
field up to polar order of = 7 and azimuthal order of 
TO = 5. Magnetic fields within the sphere are measured 
by seven linear arrays of Hall probes inserted into the 
sodium within stainless steel sheaths. These probes are 
oriented to measure either the axial or toroidal compo- 
nent of the field. Finally, two external electromagnets, 
in a Hclmoltz configuration coaxial with the impellers, 
apply a nearly uniform magnetic field throughout the 
sphere. The applied field is between and 60 G, and dom- 
inated by spherical harmonic content of £ ~ l,m = 0: the 
largest measured m ^ component of the applied field 
is less than 2% of the axisymmetric part. 



The study is conducted in the kinematic regime — the 
magnetic field is not strong enough to affect the flow. 
The strength of the Lorentz force relative to the iner- 
tial forces acting on the fluid is characterized by the 
interaction parameter (also called the Stuart number), 
N = aaB^/ pvQ, where a is the radius of the sphere, a 
and p are the conductivity and density of the fluid, re- 
spectively, and Bq and wo are characteristic magnetic and 
velocity field magnitudes. N ^ 10^^ for a total magnetic 
field of 100 G and vo = 16.0 m/s, so the magnetic field 
is not expected to alter the flow. This is confirmed by 
the linear dependence of the induced magnetic field with 
respect to the applied field. To affect the fiow we would 
expect iV w 0.1, or So ~ 180 G, a field magnitude not yet 
achieved. We note that the fluctuations, which are char- 
acterized by slower velocities, may be in a regime that is 
affected by the magnetic field. 

The axisymmetric part of the velocity field generated 
by the impellors can be expressed in cylindrical coordi- 
nates (s, 4>, z) as 

V = V$ X V<?:) + U0(s,z)^, (1) 

where $(s,2;) is the poloidal stream function. The flow 
consists of two large cells, one in the northern and one in 
the southern hemisphere. An example of this flow, based 
on measurements made in a water model of the sodium 
apparatus can be seen in the lower half of Fig. |21 
The poloidal cells flow inward at the equator and out- 
ward at the poles. The two toroidal cells flow in opposing 
directions. The flow is similar to the t2s2 flow proposed 
by Dudley and James ^ A^w which is calculated 

to magnetically self-excite at sufficiently high magnetic 
Reynolds number, Rm = fj-oaavo, where po is the vac- 
uum magnetic permeability {Rmup = pQaavup, where 
Vtip is the impeller tip speed). This study is conducted 
below the critical Rm for self-excitation, as demonstrated 
by the lack of observed growing magnetic fields. The 
Reynolds number of the fluid is Re ~ 10''; turbulent 
fluctuations of the measured flow can be as large as 20% 
of the mean, depending on location. 

Once the sphere is full of sodium the motors are started 
and a constant magnetic field is applied to the sphere. 
Hall probes sample the magnetic field at 1 kHz for 5 min- 
utes; the applied field is then subtracted from these data 
to determine the induced field. Measurements of the in- 
duced field are presented in the upper half of Fig. [3 The 
field is represented by a toroidal component, B(j,(s,z), 
and poloidal flux function, 5'(s,z), such that 

B = x\/(j) + B4s,z)4). (2) 

The toroidal magnetic field, undetectable by probes out- 
side the sphere and orthogonal to the applied poloidal 
field, is measured within the sphere by internal Hall 
probes, confirming the presence of the cj-effect. The peak 
amplitude of the toroidal magnetic field scales linearly 
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FIG. 3: Induced dipole moment versus time, for Rmup — 100 
and an applied magnetic field of GOG. IGm^ corresponds to 
13.2 G at the sphere's pole. 



with Rm, and can be larger than the magnitude of the 
applied field. 

The external induced poloidal magnetic field is decom- 
posed into its spherical harmonic components to reveal its 
spatial structure. Since the Hall probes on the sphere's 
surface lie outside regions containing currents the mag- 
netic field can be expressed as the gradient of a scalar 
magnetic potential, B = — V<i>m, which solves Laplace's 
equation. In spherical coordinates the solution to the 
potential, for the region excluding the origin, is well 
known: $™(r,0,(/)) = E^™ ^A™r-(^+i)F,"(0, ,/.), where 
Y^{6,(j)) is the spherical harmonic. The coefficients in 
the expansion, De^m, which fit the mean induced field 
are calculated using singular value decomposition. The 
induced poloidal magnetic field is predominantly axisym- 
metric; the largest components are given in Tab. The 
dominant components with £ equal to 3 and 5 are ex- 
pected due to the structure of the applied field and mean 
flow; the large measured dipole component is not ex- 
pected, as it cannot be generated by the axisymmetric 
mean flow, as will be shown below. 

Harmonic {£, m) Energy Br,max RMS Fluctuation 



1,0 (dipole) 


1.6 erg 


11.4G 


1.8 G 


2,0 


0.2 


3.1 


3.5 


3,0 


0.5 


13.6 


2.4 


4,0 


0.1 


7.1 


3.5 


5,0 


0.4 


18.3 


3.5 


1,1 


0.0 


0.8 


7.8 


2,1 


0.0 


0.6 


3.3 



TABLE I: Mean energy in the largest induced external 
poloidal harmonics, maximum mean radial field on the 
sphere's surface, and field fluctuation level for several spher- 
ical harmonic components, for Rniup = 100 and an apphed 
field of 60G. 

The induced dipole moment fluctuates dramatically in 
time around a well-defined mean, as seen in Fig. 13 Mea- 
surements indicate that the induced dipole depends on 
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FIG. 4: (a) Mean induced dipole moment versus Rmup, for 
an applied field of 60 G with a quadratic fit valid at low Rm. 
(b) Mean induced dipole moment versus applied magnetic 
field, for Rnitip — 100. A linear fit is plotted for comparison. 
Error bars are RMS fluctuation levels about the mean; the 
uncertainties in the mean values are very small (less than 
0.01 Gm'^) due to long averaging times. 



Rm (Fig.^) and upon the magnitude of the externally- 
applied field (Fig.^lD). The dipole moment's dependence 
on Rm eliminates the possibility of the measurement be- 
ing a systematic error in the analysis. The EMF depends 
linearly on the applied field, indicating that it is a kine- 
matic effect and not due to the back reaction. 

While Cowling's theorem demonstrates that self- 
excitation is not possible in axisymmetric systems, it is 
not obvious that a dipole moment cannot be induced 
by an axisymmetric velocity field exposed to an axi- 
symmetric magnetic field. The proof of this is as fol- 
lows. Consider a bounded, steady-state, axisymmetric 
system described by Eq. [3 For axisymmetric fields, 
the only non-trivial component of the dipole moment, 
/X = J X X J (fix, is oriented along the symmetry axis and 
results from currents flowing in the toroidal direction. 



SJr 



d^3 



(3) 



These currents can only be generated by the v x B force 
due to the mean fields, so using Ohm's law gives 



sJ^ = scr [v X (V^* X V(/))] • 
= a v^V^ - (v • V*) 4) 
= -crV • (v*) , 







(4) 
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where use has been made of V^" ■ </> = and the fluid 
has been assumed incompressible, V • v = 0. Inserting 
Eq. 21 into Eq. |21and making use of Gauss' theorem and 
V • n = 0, where n is the unit vector normal to the 
vessel's surface, one finds that fi^ = 0. It is interesting 
to note that it is only the dipole moment that vanishes; 
moments which include different powers of s in Eq. |2| are 
nonzero in general. This conclusion is also independent 
of geometry; any simply-connected axisymmetric system 
gives the same result. 

It is possible that an induced dipole could be gener- 
ated if mean non-axisymmetric magnetic and velocity 
field modes interacted. The stainless steel tubes which 
contain the internal Hall probes could potentially break 
the symmetry and create a mean non-axisymmetric fiow. 
However, if this were the case one would expect higher- 
order non-axisymmetric induced field components, which 
are not observed (see Tab.lJ). The mean induced dipole 
moment is present both with and without the tubes. 

Since it cannot be generated by the mean flow, the 
dipole moment must be the result of turbulence breaking 
the symmetry of the system, likely a turbulent EMF of 
some form. Any of the terms in the mean-field expan- 
sion of the EMF have the potential to yield the observed 
mean dipole moment. A toroidal a-effect could produce 
large scale toroidal currents by interacting with the ob- 
served w-effect. The small scale helicity needed for the 
a-effect might come from either a turbulent cascade or 
be produced directly by the impellers. The /3-effect leads 
to turbulent modifications of the fluid conductivity |3| . A 
nonuniform /^-effect could cause uneven distributions of 
currents to generate the dipole moment. A third possi- 
bility is the 7-effect , which expels magnetic field from 
regions of high-intensity turbulence, resulting in diamag- 
netism. The intensity of the turbulence varies with posi- 
tion, so the /3-effect and the 7-effect are both candidates 
to explain the field. 

Expanding the EMF in terms of the mean magnetic 
field may not be appropriate, since the largest fiuctu- 
ations in the magnetic field do not satisfy the scale- 
separation and homogeneity requirements usually im- 
posed in the expansion of the mean-field EMF. The 
largest turbulent magnetic fluctuations are m — 1. Their 
Gaussian probability distribution is centered at zero, con- 
sistent with a passively-advected magnetic field in a tur- 
bulent cascade of velocity fiuctuations. These m = 1 
fluctuations in B could, in principle, interact with m — 1 
fluctuations in the flow and average to give a net toroidal 
current. 

In summary, a mean dipole moment is induced in the 
experiment which cannot be produced by the mean flow. 
The induced currents are of the correct form to create 
a poloidal magnetic field, as required in the aw-dynamo 
model Q. This is the first observation of this effect in a 
laboratory experiment. Explicit characterization of the 
EMF is impossible without more detailed knowledge of 



the form of the turbulence and direct measurement of 
the fluctuating components of v and B. Future work 
will be directed towards identifying the characteristics 
of the fluctuations responsible for producing the dipole 
field. We also note that no saturation of the mechanism 
has yet been definitively observed, as might be expected 
from numerical simulations and theory [19I l2Cl| . Future 
experiments with larger magnetic fields may provide in- 
sight into the saturation mechanism. 
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